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Abstract: 

We present a next-to- leading-order QCD calculation for e + z/ e /i + z/ M jj production in 
vector-boson fusion, i.e. the scattering of two positively charged W bosons at the LHC. 
We include the complete set of electroweak leading-order diagrams for the six-particle final 
state and quantitatively assess the size of the s-channel and interference contributions in 
VBF kinematics. The calculation uses the complex-mass scheme to describe the W-boson 
resonances and is implemented into a flexible Monte Carlo generator. Using a dynamical 
scale based on the transverse momenta of the jets, the QCD corrections stay below about 
10% for all considered observables, while the residual scale dependence is at the level 
of 1%. 
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1 Introduction 



Vector-boson fusion (VBF) processes at the Large Hadron Collider (LHC) offer unique 
signatures owing to two easily identifiable forward and backward jets. This process class 
is not only useful for confirming the existence of the Standard Model Higgs boson, but 
in particular for studying its characteristics, including its couplings to both fermions and 
electroweak (EW) vector bosons [TJ[2] an d its CP properties [310]. 

VBF processes involving the scattering of vector bosons constitute an irreducible 
background to Higgs-boson production in association with two jets, in particular for 
H — > ZZ/W + W~ 41 decay modes as they share the same final states. It is therefore de- 
sirable to obtain accurate theoretical predictions and error estimates for these background 
processes. The reactions of the type pp — > VVjj — > 4Zjj + X are also seen as an important 
probe of the EW symmetry breaking itself [5]. Without the presence of the Higgs boson 
perturbative unitarity of the Standard Model at very high energy scales would be violated 
in processes involving weak- vector-boson scattering unless some other mechanism beyond 
that described by the Standard Model controls the unphysical behaviour (see e.g. Ref. [6]). 
Moreover, VBF into pairs of vector bosons is an important background to various searches 
for new physics. 

We are specifically interested in VBF processes that involve the scattering of weak 
gauge bosons and lead to final states with two jets and four leptons (charged leptons and 
neutrinos). At leading order (LO), two hard production mechanisms give rise to these 
final states. The purely EW contributions of order a 6 involve in particular the genuine 
VBF contributions, i.e. diagrams where vector bosons are emitted from the incoming 
( ant i-) quarks, then scatter and decay into pairs of leptons. In addition, there are QCD- 
production contributions of order a 4 a^ which proceed via gluon-mediated (anti-)quark 
scattering processes or processes with two external gluons and two external (anti-)quarks, 
where in both cases the two EW vector bosons are emitted from the (anti-)quark line(s). 
Since the jets in the QCD production mode tend to be closer in rapidity than in the 
EW production mode, a cut requiring a large rapidity separation between these jets or 
more generally a central jet veto suppresses the QCD production mode by two orders 
of magnitude as demonstrated in Refs. j&HH]. Owing to the different colour structure 
of QCD and EW production modes, interference terms between these mechanisms are 
doubly suppressed at LO; they only appear at sub-leading colour and if all quarks are 
identical. In our calculation, we restrict ourselves to the EW production mode. 

In this paper we focus on the process involving two jets, two positively charged leptons 
and two neutrinos in the final state, i.e. pp — > W + W + jj + X — > e+z/e/i+i^jj + X. This 
process leads to a distinct signature of same-sign high-px leptons, missing energy and 
jets. Since no gluon-initiated processes contribute to this final state at LO, it has a 
comparably low SM cross section and thus is a good candidate to search for physics 
beyond the SM. New-physics signals involving same-sign leptons originate for instance in 
i?-parity- violating SUSY models [9], in di-quark production with decay of the di-quark 
to a pair of top quarks [10J, or from the production of doubly- charged Higgs bosons [TT] . 
Moreover, it constitutes a background to double parton scattering [T2HT3] . 

Since LO cross sections carry a large uncertainty, the calculation of the NLO correc- 
tions in the strong coupling is needed to obtain a reliable prediction. For the QCD- 
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mediated contributions to pp — >■ W + W + jj + X — > e+^e/i+z^Jj + X, and later also for 
pp — > W + W~jj + X — > e+z/e^/x'jj + X, NLO results were presented in double-pole ap- 
proximation, i.e. including only diagrams with two resonant W bosons, but leptonic W de- 
cays with full spin correlations [ToT - TlT] . For W + W + the computation [15] has been subse- 
quently implemented [H] into the POWHEG BOX [T9l[20]. In a series of NLO calcula- 
tions for vector-boson scattering processes in VBF [2TH24"] , NLO results for the E W pro- 
duction mode were given for the complete process pp — > W + W + jj + X — > e+z/e/i+z/^jj + X 
in Ref. [23] including the full set of t- and w-channel diagrams (also those without res- 
onant W bosons) while neglecting s-channel diagrams and interferences between t- and 
^-channel contributions. Also this computation has recently been combined with a parton 
shower [8] using the POWHEG BOX. 

In this work we present an independent calculation of the NLO QCD corrections to 
EW W + W + jj production including leptonic W-boson decays and non-resonant diagrams. 
This constitutes the first independent check of the calculation in Ref. [23]. Moreover, we 
investigate the size of the s-channel and interference contributions at LO. 

The paper is organized as follows: In Section |2] we discuss technical aspects of the cal- 
culation, namely the organization of Feynman diagrams into building blocks (Section 12. ip 
and the evaluation of the NLO corrections (Section I2.2p . Section [3]covers numerical checks 
and comparisons with previously published results. Finally, in Section 131 we present nu- 
merical predictions for the LHC at 14TeV, while Section [5] contains our conclusions. 

2 Elements of calculation 

We have developed a general framework for the calculation of QCD corrections to vector- 
boson scattering reactions at hadron colliders [25], i.e. for EW processes of the form 
pp — > 4/jj + X with 4 arbitrary (charged and/or neutral) leptons. In this section we 
sketch the ingredients and main features of the method. 

Because the LHC experiments are conducted at TeV energies, fermion-mass effects 
are strongly suppressed and have been neglected. At the same time, only the two lighter 
generations of quarks (u, d, c and s) and leptons have been taken into account. In 
Ref. |26j, the contribution of external b quarks to Higgs production via VBF has been 
found to be at the level of 2% if VBF cuts are applied (4% without VBF cuts). For the 
processes involved in our calculation, these contributions can be expected to be of similar 
size, if not even smaller: e.g. in the W + W + channel discussed in this paper, external 
bottom (anti-)quarks would show up only accompanied by non-diagonal CKM matrix 
elements and thus be entirely negligible. Further, it can be demonstrated |23j that the 
CKM matrix can be approximated by a unit matrix provided the interferences between 
different kinematic channels as well as the s-channel contribution are negligible, which is 
verified in Section 131 

2.1 Structure of the diagrams and building blocks 

For the calculation of the LO and NLO matrix elements of the processes pp — > 4/jj + X 
we adopt a similar strategy as in Refs. [2TH24] . In order to deal with the large num- 
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(d) Type D 

Figure 1: Generic types of t-channel topologies 

ber of diagram^], we introduce generic building blocks from which the matrix elements 
can be constructed. The details of our approach differ, however, in many aspects from 
Refs. ^MM- 

For the class of processes pp — > 4/jj + X, the Feynman diagrams can be divided into 
four generic categories, taking advantage of the fact that the EW and QCD parts of 
the diagrams are largely independent of one another. Figure [T] demonstrates four generic 
types into which all i-channel diagrams involved in our calculation can be categorized. 

Type A (Figure [Taj represents the genuine VBF diagram, with two vector bosons 
radiated off the quark lines fusing in the centre to produce four leptons in the final state. 

Type B (Figure [lb]) contains two quark lines connected with a vector boson and an- 
other vector boson radiated off either of the two quark lines which decays via EW inter- 
ior uc — > dsc+^c/i+i^ there are 93 diagrams at LO. At NLO, 430 loop and 346 counterterm diagrams 
contribute, while the basic partonic real-emission process uc — > dse + v e fj, + v fl g contains 452 diagrams. 
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actions into four final-state leptons. For combinatorial reasons, 4 topologies of this type 
exist. 

Type C (Figure [Tc]) has a vector boson radiated off either one of the two quark lines 
which decays into two leptons, and two more vector bosons fusing in the central region 
to produce a second pair of leptons. Again, 4 topologies of this type exist. 

Finally, type D (Figure H~d|) sees one vector boson connecting the quark lines, and two 
more radiated either one off each quark line (diagram on the left) or both off the same 
quark line (diagram on the right) and subsequently decaying into two lepton pairs. In 
total, there are 10 different topologies of this type which are grouped together as they 
involve the same EW building blocks. 

Each generic diagram involves two QCD parts — the two quark lines with attached 
vector bosons — and one or two EW parts, namely the vector-boson scattering block and/or 
vector-boson decays into leptons. Owing to charge conservation, not all generic topologies 
give rise to Feynman diagrams once particular insertions for the external ( ant i-) quarks, 
the final-state leptons, and, correspondingly, the intermediate EW bosons are fixed. Thus, 
for instance, type B is completely absent if the final-state lepton charges add up to ±2. 

Evidently changes to either QCD or EW parts that do not alter the momenta of the 
internal vector bosons have no effect on the rest of the diagram. For instance, application 
of crossing symmetry to the upper quark line does not influence the lower quark line and 
all leptonic parts. Similar arguments hold for adding a gluon loop to either of the quark 
lines, which essentially amounts to calculating virtual NLO QCD corrections to the entire 
diagram. Since the leptonic sector of the diagram in itself can be quite complicated, it 
is advantageous to calculate these blocks only once and reuse them with different QCD 
parts. 

The calculation is performed in the 't Hooft-Feynman gauge. Diagrams with would- 
be Goldstone bosons connecting EW and QCD parts do not contribute because their 
couplings to massless fermions vanishU Factorization of parts of the diagrams can be 
achieved by inserting the polarization sums for massive vector bosons, 



={+,-,0} ''" 



for the numerators of the gauge-boson propagators coupled to a quark line, effectively 
thus cutting the diagram into blocks that can be evaluated on their own. The polarization 
vectors (k) and s* v {k) for off-shell particles are obtained by replacing the vector-boson 
mass with yW 1 in the definition of the longitudinal polarization vector, i.e. 



|k| \ 

— , cos sin 9, sin sin 6, cos 9 J . (2.2) 




The polarization vectors £±(k) do not depend on the mass and thus remain unchanged. 
Introducing 

= 4f, = --/p, (2-3) 



2 Note that would-be Goldstone bosons show up inside the EW vector-boson scattering building block 
of Figure [Tal namely in all processes involving a W + W~Z/7Z/7 vertex. 
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Figure 2: Example of a diagram split into four building blocks by applying the polar- 
ization sums (]2.ip to cut three intermediate vector bosons. 

to compactify the notation, the polarization sum ( 12.1 ft can be rewritten as 

ST = - E £ ?( k >?(k). (2.4) 

i={+,-,Q,m} 

Because of gauge invariance the contractions of fc^ with some of the building blocks vanish. 
After checking numerically that these terms do not contribute, we have omitted them in 
those type D diagrams (Figure fid)) in which two outgoing vector bosons couple to the 
same quark line since their evaluation consumes most CPU time (in comparison to the 
remaining topologies), particularly at NLO. 

Implementing the block structure by cutting all internal vector bosons that couple 
to the quark lines in the diagrams in Figure [1] allows us to not only save CPU time 
by evaluating each required block only once, but also to keep the number of required 
blocks relatively small by reusing them in multiple instances throughout all diagrams and 
even partonic processes. The diagrams we need to consider in our calculation contain up 
to three vector bosons that are being radiated off the quark lines, and the polarization 
sum has to be applied once to each of their propagators. Figure [2] illustrates how a 
diagram featuring three vector bosons can be split into four building blocks. Each splitting 
represents an insertion of one polarization sum. The resulting amplitude reads: 

- E [ A ^h,^K k Ah)] [B p 'e 3Ap ,(k 3 )} [C^'e lii) , l (k 1 )e 2 , jy (k 2 )} [D'e^M] 

i,j,k=±,0,m 



k\-M\kl-Ml 2 k\-M^ 

where l/(kf — My) are the denominator parts of the gauge-boson propagators, and 
k\ = pi — p 3 — p 5 — pe, k 2 = P2 — P4, and ks = p$ + p%, according to the notation intro- 
duced in Figure [21 

Building blocks involving leptons (shown in Figure [3]) typically involve more than one 
Feynman diagram, with the exception of the block corresponding to the EW current 
(Figure [3a]). Diagrams of type B (Figure lib"]) contain building blocks with one vector 
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(c) (d) 
Figure 3: Building blocks involving leptons. 




(c) 

Figure 4: LO building blocks involving quark lines. 

boson in the initial state and four leptons in the final state (Figure l3bl) . Building blocks 
with two external vector bosons (Figures [3c] and l3dl) are represented by a 4 x 4 array, 
each element corresponding to one term of the complete polarization sum constructed by 
cutting the two vector bosons. 

At LO, the QCD building blocks are formed by one diagram each, as shown in Figure HI 
Blocks involving two vector bosons entering diagrams of type B, C and D are repre- 
sented by a 4 x 4 array (Figure lib)) . Building blocks with three outgoing vector bosons 
(Figure l4"cl) appearing in type D are represented by 4 x 4 x 4 arrays. 

All partonic processes contributing to a specific process pp — > VVjj + X — >■ 4/jj + X 
can be constructed from up to 8 generic t-channel matrix elements listed in Table [1] 
(e.g. uc — > dse+z/g/i" 1 "^ for W + W + ) by applying crossing symmetry to reverse the flow 
of either one or both quark currents (e.g. us — >■ dce + z/ e/ u + z/^) and to construct the s- 
channel diagrams (e.g. sc — > due+z/e/i+i^), or by exchanging the outgoing lines to obtain 
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W+W+ 



uc — > dse + z/ 0/ u + z/ / 



W+W 



uc — > uce^Ueiy^n 
ds — > dse + h , e P /1 ii 
us — > use + i/ e ^/^" 
us — > dce + i/ e i> ll ii 



W"W 



ds ucz/ e e v^jj, 



W+Z: 



uc — > dce + u e fi + fi 
us — > dse + i/ e /i + yU 



ZZ: 



uc — >• uce^e - /!^ 
ds — )■ dse + e~yU + /i 
us —> use + e~fi + fi 
us — >■ dce+e'yU+yU 



W~Z: 



dc — )■ ucz/ e e yU+yU 
ds — > usz/ e e~/i + /i 



Table 1: List of generic t-channel matrix elements corresponding to the intermediate 



weak bosons produced in the fusion diagrams. Final states with less charged 
leptons are obtained by modifying the leptonic building blocks only, i.e. by 
replacing by v^v^ or e + e~fi + ^~ by v^^^v^. Moreover, W + W~ and 

ZZ in general mix if same-flavour leptonic final states like e + v e i> e e~ are con- 
sidered. All partonic processes contributing to the respective hadronic cross 
sections can be obtained from these generic matrix elements. 



li-channel diagrams (e.g. uc — > sde + z/ e /i + z/ /t ). As the order of the outgoing partons is 
obviously arbitrary, the distinction between t- and w-channel diagrams only makes sense 
if both types contribute to the same partonic process (e.g. uu — > dde + z/ e /i + z/ M ), which is 
the case if all (anti-)quarks belong to the same generation. 

With the CKM matrix approximated by a unit matrix, partonic processes which 
result from one another by interchanging all first-generation (anti-)quarks with their 
second-generation counterparts and vice versa, are described by the same matrix ele- 
ments. For instance, the partonic processes uu — > dde + z/ e /i + z/ At and cc — > sse + h , e fi + h , ^ t only 
differ in the parton distribution functions, and the matrix element can be recycled. Anal- 
ogously, partonic processes involving two different generations of quarks — for example 
sc — > ude + z/ e /i + i/ /J and du — > cse + i> e fi + i> fl — are pairwise formed by identical matrix ele- 
ments. 

In our calculation, the formulae for combining the blocks have been implemented as 
Fortran subroutines. The expressions for the individual building blocks are obtained 
by means of the FormCalc 6 package [27], which also introduces abbreviations for the 
fermion chains and thus helps to speed up the code significantly. They are evaluated 
using the Weyl-van der Waerden (WvdW) helicity formalism [28] allowing us to express 
all subamplitudes involved in the polarization sums in terms of universal WvdW spinors 
and compute them numerically. The FormCalc code is modified to transform the am- 
plitudes to the form (12. 5p . further abbreviations of spinor products are introduced, and 
each building block is exported into a Fortran module which takes the momenta and 
helicities for the particles in the building block as input. The Fortran code for each 
process is contained in a single function that can be called from within a Monte Carlo pro- 
gram and returns an array of full squared amplitudes for each partonic process, including 
all relevant colour and averaging factors. 
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2.2 Calculation of NLO cross sections and matrix elements 



Cross sections involving two initial-state hadrons at a fixed perturbative order are given 
as a convolution of the parton distribution functions (PDFs) and the partonic cross sec- 
tions dab, summed over all incoming partons resulting in contributions to the considered 
hadronic process. 

At LO, the cross section is defined by 

a pp° = 5Z / dxi / dx 2 f^°{x 1 ,fi F )ff°{x 2 ,^ F ) [ d^ m da^ b (x lPl ,x 2 p 2 ), (2.6) 

ab J0 J0 Jm 

where f a /b(%i/2, /^f) are the PDFs that give the probability to find parton a/b with a mo- 
mentum fraction xy 2 in the respective proton, da^ b (xipi, x 2 p 2 ) is the differential partonic 
Born cross section which is integrated over the m-parton phase-space $ m . While the sum 
over the incoming partons a and b is explicitly stated, summation over all outgoing parton 
configurations giving rise to non-vanishing partonic contributions to the hadronic process 
discussed is implicitly assumed. 

At NLO QCD, virtual and real corrections contribute to the cross section, which 
separately contain soft and collinear divergences. However, these infrared (IR) divergences 
cancel in the NLO cross section, if IR-safe jet observables are considered and the PDFs are 
renormalized appropriately. For mediating this cancellation, the Catani-Seymour dipole- 
subtraction technique for massless particles [2S] is applied. This procedure allows us to 
express the NLO cross section as a sum over individually finite phase-space integrals, 

<p lo =E /W NLo (xi,M F )/r°(w F ) 

„ u Jo Jo 

d$ m [dof b (xiPi, x 2 p 2 ) + da a J b {x 1 p 1 ,x 2 p 2 ) + 1 ® da^ b (x 1 p 1 ,x 2 p 2 

+ / dZi / d$ m (Kaa'(Zl) +Paa'Ol)) ® dd^^XiP^ X 2 p 2 ) 

J v 171 

+ J dz 2 J d$ m (K bb ,(z 2 ) + P bb >(z 2 )^) ®da a \,{xip 1 ,z 2 x 2 p 2 ) 
+ / d$ m+ i I da^xxpx^x^) - ^ (dVdipoie ® do- B ) ab (xipi, x 2 p 2 ) 

Jm+1 \ dipoles 

(2.7) 

with the conventions of (I2.6p . The differential partonic contributions da^ b and da\\ corre- 
spond to the virtual and real corrections, respectively. The process-independent operators 
I, K, and P are defined in Ref. [29], and <g) symbolizes the colour correlations between Born 
matrix elements and these operators (spin correlations do not appear in the given process 
class as no external gluons are involved at Born level). In the processes pp — > VVjj + X, 
the Born and virtual cross sections are built from the partonic initial states qq, qq, qq, qq, 
while in the real cross section additionally qg,gq,qg,gq contribute. 

The pro cess- dependent ingredients which are necessary for calculating the NLO cross 
sections of the processes pp — > VVjj + X are thus: 
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Figure 5: Virtual contributions Ay U and Dy corresponding to building blocks A^ v and 
in Figure [2j 

1. Born-level level matrix elements Ai B needed for d<r B , evaluated in four dimensions. 
Their construction is outlined in the previous section. 

2. One-loop virtual matrix elements Ai v needed for da v , with renormalized ultra- 
violet divergences and IR divergences regularized using dimensional regularization, 
evaluated in D dimensions. Once the QCD and EW sections of the diagrams are sep- 
arated, the transition from LO to virtual corrections can be performed by adding a 
gluon loop to either of the two quark lines. Continuing with the example in Figure EJ 
one of the building blocks A^ p or D u is replaced by A v p or D v (shown in Figure [5]) 
respectively, while the other one remains unchanged. The ultraviolet divergences 
are renormalized by adding the corresponding counterterms. Since neither of these 
changes has any influence on the overall kinematics of the diagram, the leptonic 
blocks BP and C v stay the same as in the LO. 

3. Real-radiation matrix elements Ai R needed for dcx R , evaluated in four dimensions. 
They can be created in a similar manner as in the LO case. In the example from 
Figure EJ this amounts to attaching an outgoing gluon in every possible way in the 
building blocks A w or D u and shifting the momenta ki, k<i and k^ of the intermediate 
vector bosons accordingly. The diagrams with an initial-state gluon can then be 
obtained via crossing symmetry. 

4. A set of colour projected Born-level matrix elements required to construct 
dVdipoie 8> d<j B and (K + P) <g) d<j B , evaluated in four dimensions. In the VBF ap- 
proximation of the matrix elements (no s-channel diagrams and no interferences 
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between t- and M-channel diagrams), the colour correlations turn out to be trivial 
and give rise to the same constant factor — Cp for all colour-correlated Born matrix 
elements. 

Analytical expressions for the virtual amplitudes have been generated in the same 
way as for the Born amplitudes. All divergences of the loop diagrams appear in tensor 
integrals, which are given by two, three, four and five-point functions. In our calculation, 
the tensor reduction is performed numerically in Fortran by means of the COLLIER 
library [30H35] which is based on the tensor reduction scheme developed by Denner and 
Dittmaier [3"0"|I3"T] and supports both mass and dimensional regularization scheme. 

The resonant W bosons require a proper inclusion of the finite vector-boson widths 
in the propagators. We use the complex-mass scheme, which was introduced in Ref. [36] 
for LO calculations and generalized to the one-loop level in Ref. [37]. In this approach, 
the W- and Z-boson masses as well as the Higgs-boson mass are consistently considered 
complex quantities, defined as the locations of the propagator poles in the complex plane. 
This leads to complex couplings and, in particular, a complex weak mixing angle. The 
scheme fully respects all relations that follow from gauge invariance. A brief description 
of the complex-mass scheme can also be found in Ref. [38] . 



3 Checks and comparison with existing results 

In order to verify the correctness of the calculation, comparisons with available results 
and tools have been performed at each step. The matrix elements for each partonic pro- 
cess have been verified for a set of phase-space points. In the zero-width limit, we have 
compared the Born and real-correction amplitudes for individual partonic processes with 
MadGraph 4 [39J at single precision accuracy as well as with stand-alone FormCalc 
6 [27J at double-precision accuracy and found full agreement within the numerical accu- 
racy of the calculation. For virtual corrections, the matrix elements have been generated 
with tensor reduction performed in both COLLIER and LoopTools, using mass reg- 
ularization for IR divergences. The two results agreed at the 10 -8 level. Cancellation of 
the UV divergences has been tested numerically by varying the value of the UV regulator 
£uv from 1CT 5 to 10 5 ; the resulting amplitude stays unchanged up to the level of 1CT 11 . 
Born and real matrix elements in the complex-mass scheme have been checked against 
OpenLoops [10] and found to be in full agreement within double-precision accuracy, 
both for full amplitudes (without approximations) as well as for the amplitudes in the 
so-called VBF approximation (neglecting s-channel diagrams and interferences between 
channels); in the OpenLoops framework, the VBF approximation was imposed by se- 
lecting the relevant parts of the squared matrix elements according to their particular 
colour structure. 

Furthermore, the pole structure of the virtual corrections is given by the following 
formula [21] derived from the I operator [29J: 



•M block — block 





) £m 


2 


3 " 






. fc IR 





+ C(£ir) + finite terms, (3.1) 



where A^iock and A^iock are arrays of matrix elements corresponding to the virtual 
and Born-level QCD building blocks, £ir = 2/(4 — D) stands for the IR pole, and 
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Q 2 — ~(Pi — P2) 2 = 2pi • P2 is two times the scalar product of the two quark momenta 
involved in the given building block. In the example shown in Figure [21 -Mhiock * s built 
from A^ u and Z) M , while A^^iock involves Ay 1 or D v (defined in Figure E]) instead. Relation 
(13 . 1 p has been used to verify correctness of the IR structure both for the entire virtual 
amplitude as well as for each individual building block. 

We now provide an overview of the comparisons of the full integrated cross section 
for the process pp — > e+z/c/i+z^jj + X with previously published results. All results of 
our calculation have been produced using Monte Carlo code originally developed for the 
calculation of the NLO QCD corrections to pp ->■ W + W~bb + X ->■ z/ c e + fM~ z/^bb + X [H]. 
For practical reasons — a built-in generic interface already existed — we used the tree-level 
amplitudes generated with OpenLoops [10] after cross-checking them against the ones 
obtained in the approach of Section [21 The virtual amplitudes, on the other hand, are 
constructed with the methods described in Section [21 

1. The first results for NLO QCD corrections to pp — > e + z/ e /i + z^Jj + X have been pub- 
lished in Ref. [21]. We reproduced the calculation with the same setup and in- 
put parameters. The events have been generated at the centre-of-mass energy of 

= 14TeV. In the matrix elements for all partonic processes, we neglected s- 
channel diagrams and interferences between t and u channels. The fixed-width 
scheme has been used to treat the massive propagators, with the exception of the 
Higgs couplings where M z and M w have been kept complex due to technical reasons, 
while sin 8 W and cos 6 W are real. The factorization and renormalization scales have 
been set to ji-p = /U R = M w . The values of the VBF cuts are taken from Ref. [24J, 
however, the requirement that the charged leptons fall between the tagging jets in 
rapidity (I4.15P has been omittedjf] 

The LO and NLO results for the two PDF sets used in Ref. [21] are shown in Table [21 
For both LO and NLO cross sections, the relative deviation between the results of the 
two calculations is only ~ 0.2% or even smaller. These small discrepancies could 
be attributed to the slight differences in applying the width scheme (see above). 
However, assuming a statistical error of the results of Ref. [21] of per-mille (it is 
stated to be at the sub-per-mille level), the difference amounts to only 2a and is 
thus acceptable. The differences between the two PDF sets are at the level of 5% 
at LO and of 2% at NLO. 

2. In Ref. [8] the results for pp — > e + i/ c /i + z/^jj + X have been presented at the centre- 
of-mass energy of \fs = 7 TeV. While the main focus of Ref. [8] lies on the inclusion 
of parton-shower effects, the NLO QCD result for the cross section is also shown. 

As in the previous case, we have reproduced the computation of the cross section 
with the same setup, parameters, and kinematic cuts. The factorization and renor- 
malization scales have been set to a dynamic value defined as 



/iR = ^ F = PTJl+PTj2+ 5 T ' Wl+ET ' W2 with E T , Wl = . fj^Tti 



/ 2 - y iyi w t^t,w 1/2 - 

(3.2) 



3 Private communication during comparisons revealed that in the results presented in Ref. [24j this cut 
has been omitted. A corrected version of the article can be found at arXiv:hep-ph/0907.0580 
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PDF set cr[fb] 


cjoz[fb] 


5[%] 


Leading order 


CTEQ6L1 1.4746(7) 
MSTW08 1.4061(7) 


1.478 
1.409 


-0.23(5) 
-0.21(5) 


Next-to-leading order 


CTEQ6M 1.405(1) 
MSTW08 1.372(1) 


1.404 
1.372 


+0.10(9) 
-0.00(9) 



Table 2: Comparison of the integrated cross section a with the cross section ctjoz 
presented in Ref. |24| for the W + W + production processes at NLO. The 
error estimates for a are shown in brackets and affect the last digit of the 
result. The statistical error of the cross section ojoz is stated to be at the 
sub-per-mille level and is not taken into account in the last column. 



Here, Pt,w 1/2 represents the transverse momentum of the respective same-flavour 
lepton-neutrino pair, and pxj; are the transverse momenta of the two tagging jets. 
This choice of scale is slightly different from the one in Ref. [5] where, as required by 
POWHEG, the jets of the underlying Born process were used. Another difference 
between the two calculations lies in the width scheme; while our calculation used 
the complex-mass scheme, the results in Ref. [8] have been obtained using the fixed- 
width scheme. The impact of the different scheme choices, however, is known not 
to exceed the level of a few per-milie here. 

The results for the total NLO cross section shown in Table |3]for the two calculations 
differ by 2.3%. Considering the small differences in the scale choice, and in particular 
the statistical error of Ojz, the level of agreement is fully acceptable. 



PDF set <r[fb] 


<7 JZ [fb] 


5[%] 


Leading order 


MSTW08 0.16836(8) 


/ 


/ 


Next-to-leading order 


MSTW08 0.1961(2) 


0.201(3) 


-2.3(1.5) 



Table 3: Comparison of the integrated cross section a with the cross section ajz pre- 
sented in Ref. [8J for the W + W + production processes NLO. The error esti- 
mates are shown in brackets and affect the last digit (s) of the result. 
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4 Numerical results 



4.1 Input parameters and setup 

All EW Standard Model parameters used in the calculation are determined from the 
values of the Z-boson mass Mz, the W-boson mass Mw, the Higgs-boson mass Mh, and 
the Fermi coupling constant [32] ■ The EW mixing angle 6 W is defined as 



cos 6L 



M 



w 



My 



(4.1) 



The fine-structure constant a is evaluated from Gn, Mw and Mz according to 



a 



V2M^G^ 



71 



1 - 



M 



w 



Ml 



(4.2) 



which takes into account dominant effects associated with the running of a from zero to 
the W-boson mass and absorbs leading universal corrections oc G^m^ associated with the 
p parameter [33]. 

For all results presented in this section, we make use of the PDF set MSTW2008 [33] , 
i.e. MSTW2008LO and MSTW2008NLO for LO and NLO cross sections, respectively. 
Throughout, the NLO value of the strong coupling constant a s provided by this PDF set 
is used (no strong couplings appear at the LO). 

The decay widths of the unstable intermediate vector bosons are calculated at NLO 
QCD level according to 




where I runs over all charged leptons and neutrinos, q runs over the five light quarks, 
N c = 3 is the number of quark colours, and Qi, Q q and If, I? are the charges and third 
isospin components of the respective leptons and quarks. As the leptonic decays of the 
EW bosons do not receive QCD corrections at NLO, we may use the same NLO values of 
the widths, provided by (14.31) . both at LO and NLO without introducing inconsistencies 
in W or Z branching ratios. 
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Throughout the subsequent numerical discussion, we evaluate cross sections and dis- 
tributions to pp — > e + z/ e /i + ^jj + X at the centre-of-mass energy a/s = 14 TeV, using the 
following Standard Model parameters, 

M w = 80.399 GeV, T w = 2.099736097449861 GeV, 

M z = 91.1876 GeV, T z = 2.509659634331562 GeV, 

M H = 125 GeV, T H = 4.07 x 10" 3 GeV 

G F = 1.16637 x 10~ 5 GeV -2 , a s (M w ) = 0.1225519862138941, (4.4) 

where a s (M w ), T w , and T z are calculated values and thus stated at machine precision to 
facilitate comparisons with our results. The decay width of the Higgs boson Th depends 
on the chosen mass, and its value is taken from Ref. [45J. 

To treat the propagators of the unstable massive intermediate particles (W, Z, and 
Higgs boson), we use the complex-mass scheme [36H38] . in which the masses are globally 
replaced according to 

rCMS _ . [77 



My = \jM y - iMyTy. (4.5) 

Complex masses are then introduced everywhere in the Feynman rules, including the weak 
mixing angle, 

2 a v ~ iM w r w 

cos e - "> Ml - iM z r z ' (46) 

rendering the couplings complex. Note, that real masses and mixing angle are used to 
determine the input values ( 14. 2 p and (14. 3p . 

The LO cross section has been evaluated for three different setups. In the first, we 
only take into account the t-channel and w-channel diagrams and completely disregard 
the interferences between them. This approximation corresponds to the setup in Ref. [21] 
and is referred to as VBF approximation. In the second setup we include interferences 
between t-channel and it-channel diagrams, and in the third one we calculate the complete 
cross section including t-, u- and s-channel diagrams and all interferences. This allows to 
assess the size of the s-channel and interference contributions. For the NLO cross section, 
the s-channel diagrams and interferences are neglected throughout, both in virtual and 
real corrections. 

As in the comparisons of integrated results in the previous section, all cross sections 
and distributions have been produced using Monte Carlo code developed for the cal- 
culation of the NLO QCD corrections to pp ->■ W + W~bb + X -> u e e + ^"z/^bb + X [4"T] . 
using adapted tree-level amplitudes generated with OpenLoops [ID] while the virtual 
corrections were calculated according to the method described in Section [2J 



4.2 Jet recombination and phase-space cuts 

In order to enhance regions of the phase-space where VBF-type processes can be observed 
experimentally and QCD background is reduced, a number of kinematic cuts has been 
imposed at the Monte Carlo level. These cuts are selected to focus on the phase-space 
region dominated by VBF processes, which typically contain two hard jets with large 
rapidity separation, the so-called tagging jets, and most of the decay products of the 
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vector bosons in the central detector region. Further, a set of lepton cuts is applied to 
ensure that the charged leptons, which define the respective final state, are well-observable 
and separated from the jet activity. The set of cuts to be precisely defined in this section 
follows the proposal of Ref . [21] . 

To be considered protojets which eventually give rise to hadronic jets in the final state, 
outgoing QCD partons have to fulfil the requirement 



1, P0+Pz 

2 Po-Pz 



< 5. 



(4.7) 



Starting with these protojets, the jet reconstruction is performed using the &t algo- 
rithm [1HJI17] with the resolution parameter D — 0.7. In order to be clearly distinguished 
from QCD background, the resulting jets must satisfy the transverse momentum and 
rapidity cuts 



pi ,+pl .>20GeV, \ m \ 



1 ln Poj + Pzj 

2 Poj - P 



< 4.5. 



(4.f 



At least two jets have to pass this criterion, and the two jets with the highest transverse 
momenta are denoted as tagging jets, on which the following additional restrictions are 
imposed. The two tagging jets must have a minimum invariant mass, 



Mn = V(P0Ji +Poj 2 ) 2 - (Ph +Pj 2 ) 2 > 600 GeV, 
they must be located in the opposite hemispheres of the detector, 

Vh x Vh < °> 

and show a large rapidity separation, 

A %j = \Vh-Vh\ > 4 > 



(4.9) 



(4.10) 



(4.11) 



in order to further suppress the gluon-induced production mode and background processes. 
The charged leptons are required to pass transverse-momentum and rapidity cuts, 



p T ,i>20GeV, |ja|<2.5. 



(4.12) 



To ensure that they are well separated from one another and from the two tagging jets, 
we impose the additional cuts 

ABji > 0.4, ARa > 0.1, (4.13) 

where the quantity ARij is a measure of distance in rapidity and azimuthal angle, defined 
as 

(4.14) 



ij 



yj (Vi ~ Vj) 2 + S - <t>jY 



where yi, yj and fa, <pj are the rapidities and azimuthal angles of the respective particles. 
Finally, the rapidities of the charged leptons are required to fall between the tagging-jet 
rapidities, 

2/j min < V\ < 2/jmax, (4-15) 

which again points out a typical feature of the VBF-production mode. 
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Figure 6: Scale dependence of the LO (dotted blue line) and NLO (solid red line) cross 
section for the fixed (Figure [6a]) and dynamic scale (Figure I6t>|) as a function 
of the scale parameter £. 



4.3 Integrated cross sections 

We have chosen two types of scales to demonstrate their effects on the behaviour 
of the NLO distributions for selected observable quantities. In the fixed-scale (FS) 
choice, both factorization and renormalization scales have been set to the mass of 
the W boson, which sets a natural scale for the total cross section of the process 
pp — > W + W + jj + X — > e + z/ e /i + z/ /i jj + X, and varied by a factor £ around this central value, 

Hf = fiR = £M W - (4.16) 

Since this FS choice turns out to result in strongly phase-space dependent K factors — in 
particular in the high-energy tails of distributions — a dynamical scale (DS), 

fiF = m = i^/pvji -Ptj 2 - (4.17) 

has been considered as well. Unlike the DS chosen in Ref. [21], (14. 17ft only depends on 
final-state momenta and can thus be easily defined in an IR-safe way also at NLO. The 
scale in (I4.17P has been chosen to flatten the variation of the K factor in the high-energy 
tails of pr,j as well as of other energy- dependent distributions, which is demonstrated in 
the next section. 

The dependence of the total cross section on the parameter £ for both scale choices is 
depicted in Figure E] for a variation of £ in the range 1/8 < £ < 8. In the conventionally 
chosen range 1/2 < £ < 2, the scale variation of the LO cross section which only depends 
on £ via /xp (as /Jr, only enters in cv s (//r)) is about ±10%, while at NLO it is reduced to 
about ±2% of the total cross section for the FS choice and ±1% for the DS choice. For 
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IU11 L J 


rr LO f fhl 

v i y r -pinx l j 


^VB F [fb] 


^VBF[fb] 

V ur L J 


1/8 


1.6763(2) 


1.6755(2) 


1.6771(2) 


1.198(2) 


1/4 


1.4956(2) 


1.4949(2) 


1.4964(2) 


1.264(1) 


1/2 


1.3467(2) 


1.3461(2) 


1.3474(2) 


1.2903(9) 


1 


1.2224(2) 


1.2218(2) 


1.2230(2) 


1.2917(8) 


2 


1.1173(2) 


1.1168(2) 


1.1179(2) 


1.2778(7) 


4 


1.0275(2) 


1.0270(2) 


1.0280(2) 


1.2544(6) 


8 


0.9499(2) 


0.9494(2) 


0.9504(2) 


1.2253(6) 



Table 4: Integrated cross sections for LO including all channels and interferences 
(°fuu)' ^ or including t-u interferences but neglecting s-channel diagrams 
(^"VBF+int)' ^ or LO in the VBF approximation, i.e. neglecting all s-channel 
diagrams and interferences (ovbf)> and for NLO in the VBF approximation 
(°vbf)- ^he integration-error estimates are shown in brackets and affect the 
last digit of the respective result. 



scales down to £ = 1/8, the scale dependence for the DS is less pronounced (~ —7% of 
the total cross section at £ = 1) than in the case of the FS (~ —20% of the total cross 
section at £ = 1). For both scale choices the maximum of the NLO curve is located in 
the vicinity of £ = 1 (reflecting a small residual scale dependence in this region), so this 
value is chosen for the subsequent numerical discussions. For £ = 1, the overall K factor, 
denned as a NLO /a LO , is 

K FS = 0.976, K DS = 1.056, (4.18) 
respectively, for the two scale choices under consideration. 

The dedicated VBF cuts listed in Section 14.21 prefer t- and w-channel kinematics, 
whereas s-channel configurations are strongly suppressed by the requirement of final-state 
jets with large rapidity separation and invariant mass. Moreover, interferences between t 
and u channels, showing up in partonic processes with identical final-state ( ant i-) quarks, 
are suppressed by the condition that the tagging jets have to be located in opposite- 
forward and backward — regions of the detector. It can therefore be argued [24J that the 
s-channel and interference contributions can safely be neglected if the VBF cuts are ap- 
plied. In order to verify this claim, the LO cross section has been evaluated for three 
different sets of matrix elements; the results (obtained using the DS with selected val- 
ues of parameter £) can be found in Table HJ Here, cx^ stands for the cross section 
that includes all channels and interferences, while c" V BF+int contains the complete t- and 
it-channel contributions with interferences but no s-channel contributions, and <Jy^ F con- 
tains only squares of t- and w-channel contributions but no interferences. One can see 
that for the cross section within our set of VBF cuts the effect of the t-u interferences is 
at the level of —0.05%, and the contribution of s-channel diagrams at the level of +0.1%. 
This confirms that cryg F can be considered a very good approximation of the full LO 
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cross section. For this reason, the NLO cross section has been evaluated using only t- 
and li-channel contributions without interferences between them in order to improve the 
speed of the calculation. The values of OyBF for different values of £ can be found in the 
fifth column in Table HJ 

4.4 Jet distributions 

All distributions shown in this and the following section are evaluated in the numerical 
setup of Section 14.14 using the kinematic cuts introduced in Section I4.2( the scale choice 
applied is stated in each case. 

Two plots are presented for each observable: the one on the left depicts the LO 
and NLO predictions, the uncertainty of which is indicated by error bands resulting 
from variation of the given scale within 1/2 < £ < 2, while the plot on the right shows 
the LO and NLO predictions normalized to the LO result at the central scale, i.e. 
Klo(Q = da L o(0/ dor Lo(£ = 1) (dotted blue line), andX N Lo(0 = dcr N Lo(0/ dor Lo(£ = 1) 
(solid red line). The blue band in this case corresponds to the relative scale uncertainty of 
the cross section at LO, and the central curve of the red band represents the conventional 
K factor K NLO (£ = 1). 

Figure [7J shows the LO and NLO cross sections as functions of the transverse momenta 
of the harder (in terms of pj) of the two tagging jets in the range px,j max < 400 GeV. 
Figure [7a] displays the dependence for the fixed scale (I4.16P and Figure [7b] for the dy- 
namic scale (14.171) . In both cases, the distribution peaks at px ~ HO GeV, confirming the 
preference of the high-px regions by the VBF tagging jets, while the probability to find 
a jet at lower values of px is slightly larger at NLO than at LO. One can observe that 
^(Prjmax) grows noticeably in low-px regions for both FS and DS towards the value 1.3, 
while in the larger-px regions it drops to 0.8 in case of the FS (Figure ITaj) and remains 
very close to 1 for the DS (Figure [7b]) . which is a behaviour that motivated the choice of 
the DS in the first place. 

A similar behaviour can be observed with the transverse momentum of the softer tag- 
ging jet Pxj min , as shown in Figure [8] Here, the peak of the distribution is at px ~ 60 GeV, 
indicating that the cut of 20 GeV on the transverse momenta of the tagging jets does not 
impose any significant reduction to the overall cross section. The variation of i^(px,j min ) 
is less pronounced than in the case of Px,j max > while the choice of DS again shows an 
improvement at reducing contributions from higher-order corrections. 

The rapidity of the tagging jets is another distinguishing feature of the VBF processes, 
as they exhibit very little jet activity in the central region. Absolute rapidity distributions 
for the harder and softer (in terms of px) tagging jets are shown in Figures 191 and [TUl 
respectively. One can see that the probability to find the harder jet peaks at absolute 
rapidity of y ~ 2.6 while the softer jet is most likely to be found with absolute rapidity 
of y ~ 3.1. This is in sharp contrast to the behaviour for the QCD production mode for 
W + W + , where the jet rapidity peaks at 0, dominating the central rapidity region [8]. This 
production mode thus can be suppressed dramatically by imposing a cut on the separation 
of individual jet rapidities Ayjj (14. lip . One can see from Figures 191 and [TOl that for both 
scale choices the rapidity-dependent K factor for the hard jet ^(|yj hard |) has a tendency 
to grow for large values of rapidity. This might be attributed to the fact that while only 
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Figure 7: Transverse momentum distribution of the tagging jet with the higher 

for the fixed (Figure I7ap and dynamic scale (Figure I7bp on the left and the 
corresponding K factor represented by the solid (red) line on the right. 



two final-state partons are present at LO, in NLO the tagging jets are selected from up 
to three partons which might lead to greater dispersion in the rapidity distribution (see 
discussion in Ref. [21] )• As in the case of the transverse-momentum distributions, the DS 
shows an improvement over the FS in the variation of the K factor (Figures [9b] and UObD . 

At hadron colliders, QCD processes typically occur at smaller energy scales than EW 
processes. Due to the back-to-back geometry and large momenta of the tagging jets in 
VBF, the invariant mass My denned in ( 14. 9 p can easily exceed 1 TeV, which is typically 
not the case for any QCD background process, particularly if incoming gluons, which 
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Figure 8: Transverse momentum distribution of the tagging jet with the lower pr for 
the fixed (Figure I8a|) and dynamic scale (Figure l8bj) on the left and the 
corresponding K factor represented by the solid (red) line on the right. 



prefer smaller momentum fractions than valence quarks, are involved. For this reason, 
invariant-mass cuts are applied to distinguish these types of processes. Figure [TT1 provides 
the distribution for the invariant mass of the tagging jets for both FS and DS. The peak 
of the distribution is located at approximately 1100 GeV, both at LO and NLO. On a 
qualitative level, the behaviour of the distributions as well as the K factor K(M jJ ) is in 
good correspondence to those shown in Ref. [21], in particular for the FS which is set to 
Mw in both cases. While the DS in Ref. [24J is set to the momentum transfer between 
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Figure 9: Absolute rapidity distribution for the harder tagging jet for the fixed 
(Figure l9a|) and dynamic scale (Figure l9b|) on the left and the corresponding 
K factor represented by the solid (red) line on the right. 

the incoming and outgoing partons rather than to a /Ptj7 t Ptj7; it nas a similar effect on 
the behaviour of the NLO distribution, reducing the scale variation of the K factor. 



4.5 Leptonic distributions 

The decay products of the intermediate gauge bosons in VBF processes can be found 
almost exclusively in between the tagging jets, in the central region of the detector. 
This kinematic feature is used to further suppress in particular irreducible background 



21 




pp -> fee+f^AJ+jj+X @ = 14TeV 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 

W'soft I 

(a) fi F = fj, R = M w 




<Wc% soft |[fb] 
0.7 




K 
1.3 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 

(b) /i F = HR = y/PTj! -PT,32 



Figure 10: Absolute rapidity distribution for the softer tagging jet |yj soft | for the fixed 
(Figure flOap and dynamic scale (Figure [l0b|) on the left and the corre- 
sponding K factor represented by the solid (red) line on the right. 



from gluon-mediated contributions to the process pp — > W + W + jj + X — > e + i/ e /i + z/ M jj + X, 
which do not show these characteristics, e.g. by imposing lepton-rapidity cuts like the one 
in (Q5]l . 

The high leptonic activity in the central region can be well observed in the rapidity 
distribution of the harder (in terms of pt) charged lepton, shown in Figure [12j The 
distribution shows a preference for the rapidities close to zero, while decreasing quickly as 
the values approach those of the tagging jets. One can see that for both the FS and the DS 
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Figure 11: Distribution in the invariant mass of the two tagging jets Mjj for the fixed 
(Figure [l lap and dynamic scale (Figure fllbp on the left and the corre- 
sponding K factor represented by the solid (red) line on the right. 

the K factor remains quite flat. As a similar behaviour can be observed in all presented 
leptonic distributions, the following distributions are only shown for the DS (I4.17p . 

The distribution of the transverse momentum for the harder charged lepton and of the 
missing corresponding to the vectorial sum of the transverse momenta of the electron 
neutrino and the muon neutrino from the W decays are shown in Figures [TBI and fT^t 
respectively. The K factors decrease with increasing transverse momentum and are close 
to 1 for large pr- 
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Figure 12: Rapidity distribution for the harder charged lepton for the fixed 
(Figure f!2a|) and dynamic scale (Figure [l2bp on the left and the corre- 
sponding K factor represented by the solid (red) line on the right. 



Angular distributions of the decay products of the vector-boson pairs produced in VBF 
processes are of particular interest to the SM Higgs searches at colliders as the leptons 
have a tendency to fly in the same direction in case of a Higgs signal jHlll9] . This is due 
to the fact that a scalar particle decays into a pair of W + W~ which subsequently decays 
into charged leptons and neutrinos, with the leptons preferably close to each other due to 
the left-handed structure of the EW force. This is not the case in the W + W + production 
processes where a SM Higgs decay into both intermediate vector bosons in question is pro- 
hibited by the charges. However, a similar situation could arise in the presence of a doubly- 
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Figure 13: Distribution of the transverse momentum of the harder charged lepton for 
the dynamic scale on the left and the corresponding K factor represented 
by the solid (red) line on the right. 
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Figure 14: Distribution of the missing transverse momentum produced by two outgo- 
ing neutrinos for the dynamic scale on the left and the corresponding K 
factor represented by the solid (red) line on the right. 



charged scalar resonance, which would be produced in the VBF mode and subsequently de- 
cay into W + W + — >■ e + z/ e /i + i/ M . The process pp — > W + W + jj + X — > e + h , e fi + h , fM ]} + X would 
then deliver the dominating irreducible background, and the angular distributions of the 
decay products should be well known to allow for an efficient background suppression. 
As demonstrated in Figure [TB] where the azimuthal angle <f) e + IJi + separating the charged 
leptons e + and fi + in the plane transverse to the beam direction is depicted, the final-state 
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Figure 15: Distribution of the azimuthal angle between the charged leptons e + and 
/i + for the dynamic scale on the left and the corresponding K factor rep- 
resented by the solid (red) line on the right. 



leptons are located preferentially in opposite directions in the azimuthal plane. Similarly 
as with other leptonic observables, the NLO corrections have only a modest effect in this 
distribution, as exhibited by a flat K factor. 



5 Conclusions 

This article presents a method for evaluating the NLO QCD corrections to the electroweak 
production mode of processes of the form pp — >■ 4Zjj + X, associated with vector-boson fu- 
sion of an intermediate vector-boson pair, including both resonant contributions where the 
final-state leptons are produced via vector-boson decay as well as non-resonant ones. The 
Feynman diagrams are divided into independent building blocks using internal polariza- 
tion sums and evaluated with the FeynArts + FormCalc package in Mathematica 
using the Weyl-van der Waerden helicity formalism. The block structure separates the 
electroweak and QCD sectors of the diagrams, allowing one to apply the QCD corrections 
only to building blocks involving quark lines while electroweak building blocks are eval- 
uated merely at tree level, improving thus speed of the resulting Fortran code. The 
phase-space integration is performed by a multichannel Monte Carlo generator imple- 
mented in CH — h which allows to calculate arbitrary distributions. 

The described method is applied to the process pp — > e+z/e/i+i/^jj + X, which is inter- 
esting in its own right as well as constitutes a background to many collider searches both 
within and beyond the Standard Model. The numerical analysis is performed using typ- 
ical vector-boson-fusion cuts chosen to enhance contributions of the vector-boson-fusion 
kinematics and to suppress QCD background. The impact of the s-channel diagrams and 
interferences between t and u channels is analyzed at LO and found to be entirely negligi- 
ble within vector-boson-fusion cuts. The NLO corrections turn out to be around 5% of the 
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LO cross section. The renormalization- and factorization-scale dependence of the cross 
section at NLO is reduced, amounting to about ±2% for the fixed scale \i = £M W an d 
±1% for the dynamical scale \i = £a/Ptji • PT,j 2 > w h en varying £ within 1/2 and 2 about 
£ = 1, while the LO results change significantly, by ±10%, in the same range. Further- 
more, a set of kinematical distributions for jets and final-state leptons has been presented, 
demonstrating the effects of the NLO corrections and the impact of the two scale choices. 
While the fixed scale results in a strongly decreasing K factor in the high-energy tails 
of the distributions, the K factor for the dynamical scale approaches one in these kine- 
matical regions. Using the dynamical scale at leading order provides an approximation 
to the next-to-leading-order result with an accuracy below about 10% for all considered 
distributions. 
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